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We study nonclassical correlations beyond entanglement in a family of two-mode non-Gaussian 
states which represent the continuous-variable counterpart of two-qubit Werner states. We evaluate 

quantum discord and other quantumncss measures obtaining exact analytical results in special 
instances, and upper and lower bounds in the general case. Non-Gaussian measurements such as 
photon counting arc in general necessary to solve the optimization in the definition of quantum 
discord, whereas Gaussian measurements are strictly suboptimal for the considered states. The gap 
between Gaussian and optimal non-Gaussian conditional entropy is found to be proportional to a 
measure of non-Gaussianity in the regime of low squeezing, for a subclass of continuous-variable 
Werner states. We further study an example of a non-Gaussian state which is positive under partial 
transposition, and whose nonclassical correlations stay finite and small even for infinite squeezing. 
Our results pave the way to a systematic exploration of the interplay between nonclassicality and 
non-Gaussianity in continuous-variable systems, in order to gain a deeper understanding of — and 
to draw a bigger advantage from — these two important resources for quantum technology. 



PACS numbers: 03.67.-a, 03.65.Ta, 42.50.Dv 
I. INTRODUCTION 



Two decades after pioneering contributions of the likes 
of quantum cryptography [1] and teleportation [2] , quan- 
tum information science has nowadays acquired a status 
of maturity [3] , having witnessed theoretical and experi- 
mental advances which demonstrated in several ways the 
power of quantum technology [4]. Nonetheless, funda- 
mental questions of inherent practical relevance remain 
open, broadly concerning the proper identification of the 
ultimate resources behind such a power. Paradigmatic 
instances of systems suitable for quantum protocols have 
been, for instance, registers of pure qubits for discrete- 
variable ("digital") quantum computation and cryptog- 
raphy [3, 5], and multimode Gaussian states of radia- 
tion fields for continuous- variable ("analog") quantum 
communication and information processing [6]. Entan- 
glement [7] has been the crucial resource on which most 
studies have focused so far as its presence is crucial for 
better-than-classical communication performances. Re- 
cently, however, there has been more of a push for re- 
search in quantum information to go beyond the bound- 
aries of its first generation. In particular, the need to 
consider more realistic setups where mixcdness affects 
quantum computations in the discrete variable frame- 
work [8] , and the need to go beyond the nutshell of Gaus- 
sian states and operations to achieve universality in con- 
tinuous variable computation [9] have both been strong 
motivations. Nonclassicality (i.e. the quantumncss of 
correlations beyond and without entanglement) [10] and 
non-Gaussianity [11] are currently under the limelight 
due to being recognized as "power-ups" for quantum 
technology, in particular in applications such as quantum 



computation [8, 9], quantum communication [12, 13] and 
metrology [14, 15]. 

In this paper we investigate the nonclassicality of cor- 
relations in continuous-variable (CV) systems beyond 
Gaussian states. Nonclassical correlations, including and 
beyond entanglement, can be quantified for instance by 
the quantum discord [16, 17], a measure that aims at cap- 
turing more general signatures of quantumncss in com- 
posite systems. They are associated e.g. with the non- 
commutativity of quantum observables and with the fact 
that local measurements generally induce some distur- 
bance on quantum states, apart from very special cases 
in which those states admit a fully classical description. 
The interest in quantum discord has risen significantly 
since it was recognized as the potential resource behind 
the quantum speed-up in certain mixed-state models of 
quantum computing such as the DQCl [8], in which en- 
tanglement is negligible or strictly vanishing [18]. Quan- 
tum discord [16], along with other similar nonclassicality 
measures [19] such as the (ameliorated) measurement- 
induced disturbance [20-22] and distance-based quanti- 
fiers such as the geometric discord [23] and the relative 
entropy of quantumncss [24, 25], typically feature; non- 
trivial optimizations in their definitions. In particular, it 
is necessary to identify the least disturbing measurement 
to be applied on one or more subsystems to extract those 
nonclassical correlations, rendering their exact computa- 
tion a formidable task. Closed formulae are available for 
the quantum discord of a special subclass of two-qubit 
states [26] and of general two-mode Gaussian states un- 
der the restriction of Gaussian measurements [27]. 

We consider here a family of two-mode states that 
are the CV counterparts of two-qubit Werner states 
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[28]. They are non-Gaussian states obtained as mix- 
tures of two Gaussian states [29] , an entangled two-mode 
squeezed state and a two-mode thermal product state, 
and find applications in C V quantum cryptography [30] . 
Studying their nonclassicality beyond entanglement is 
particularly interesting from a fundamental point of view, 
as they offer a unique testbed to compare the role of 
Gaussian versus non-Gaussian measurements to extract 
correlations with minimum disturbance from general two- 
mode CV states. Gaussian states and operations are 
known to satisfy sharp extremality properties [31] in the 
space of all CV states. However, such results do not 
apply to quantum discord, rendering the understanding 
of the structure of (non) classical correlations even sub- 
tler in CV systems compared to finite-dimensional sys- 
tems, for which general results are known instead [32]. 
Wc show that the states analyzed here represent in- 
stances of quantum states carrying genuine CV nonclassi- 
cal non-Gaussian correlations, for which the optimization 
in the discord is achieved only by an infinite-dimensional 
component, hence genuinely CV, non-Gaussian measure- 
ment.^ In the following we prove that for a particular 
type of the non-Gaussian CV Werner state [29] (obtained 
by mixing a two-mode squeezed state with the vacuum) 
quantum discord can be computed analytically. The 
global optimization in its definition is achieved exactly 
by photon counting, and such a measurement is strictly 
less disturbing than any Gaussian measurement. The gap 
between optimal Gaussian measurements (homodyne de- 
tection) and optimal non-Gaussian measurements (pho- 
ton counting) is quantified by the difference of the asso- 
ciated conditional entropies. The gap is then compared 
with an entropic measure for the non-Gaussianity of the 
states [33], and found to be proportional to it in the low 
squeezing regime. The exactly computed discord turns 
out to be equal to other measures for the quantumness of 
correlations as well in this special case [22, 25]. Next, we 
analyze discord in the most general CV Werner state, ide- 
ally complementing the analysis of the two-qubit Werner 
state performed originally in [16] (see also [26]). We show 
that photon counting provides in general a (not necessar- 
ily optimal) upper boimd on quantum discord, that co- 
incides with the measurement-induced disturbance [20], 
and we derive a lower bound on discord as well. Finally, 
we give an example of a CV state which is positive un- 
der partial transposition (PPT [34]), and show that it 
carries "weak" nonclassical correlations, signaled by an- 
alytically computable upper and lower bounds on discord 
which are close to each other and stay small and finite 



^ One might argue, otherwise, tliat e.g. any two-qubit state is an 
example of a non-Gaussian state; for them, however, the compu- 
tation of diseord involves two-component measurements, which 
are certainly non-Gaussian but span just a two-dimensional 
space. Here we focus instead on CV states and measurements 
spanning the whole infinite-dimensional Hilbert spaces, and not 
just finite-dimensional truncations of it. 



even for infinite squeezing. Our results provide insights 
into the relationship between the general quantumness of 
correlations, entanglement distillability and separability 
in CV systems outside the Gaussian scenario. From a 
practical perspective, our results identify the key role of 
non-Gaussian measurements such as photon counting to 
access and extract all nonclassical correlations in general 
CV states, even in the particular case of non-Gaussian 
states with a positivc-everywherc Wigner function such 
as the mixtures of Gaussian states studied here. 

The paper is organized as follows. In Sec. II we set up 
notation and recall the definitions of quantum discord 
[16] and of the CV Werner states [29]. In Sec. Ill we 
study a special type of CV Werner state and calculate 
its discord (and other nonclassicality indicators) exactly, 
proving that photon counting constitutes a globally opti- 
mal measurement strategy, and issuing comparisons be- 
tween the degree of non-Gaussianity [33] and the amount 
by which Gaussian measurements fail to be optimal on 
the considered class of states. In Sec. IV we address the 
general case and provide exact lower and upper bounds 
for the quantum discord of arbitrary CV Werner states. 
In Sec. V we study a special type of CV state obtained 
by partial transposition from a CV Werner state, and 
investigate its degree of nonclassical correlations in dif- 
ferent entanglement regimes. We draw our conclusions 
in Sec. VI. 



II. BASIC NOTATION AND DEFINITIONS 
A. Quantum discord 

Quantum discord is a measure of the nonclassicality of 
correlations in a bipartite quantum state pab defined as 
a difference [16] 

T^ipAs) = IqipAs) - J{pab) (1) 

between the quantum generalizations of two classi- 
cally equivalent expressions for the mutual information, 
namely the quantum mutual information 

Xq(pAB) = S{pa)+S{pb) - S{pab), (2) 

where S{p) = — Tr(plnp) is the von Neumann entropy, 
and the so called one-way classical correlation [17] 

J{Pab) = S{pa) - mf U{u,}{A\B), (3) 

with 'H{Yi.}{A\B)=Y,iPB{i)S{pA\i) being the condi- 
tional entropy of A given a positive operator valued mea- 
surement (POVM) {WB{i)} fias been performed on B. 
Here pa.b = Tr^.^ [pab] are reduced states of subsystems 
A and S, respectively, pA\i = ^'i:B\^B(i)pAB]/PB{i) is 
the conditional state obtained upon detecting the POVM 
element Yi.B{i) on B sndpsii) = T^\^-B{i)pAB\ is the cor- 
responding probability. Without loss of generality one 
can restrict to rank 1 POVMs in what follows [35]. 
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Subtracting Eq. (3) from Eq. (2) one gets for quantum 
discord the following formula: 



S{pab)+ inf^n[n,}{A\B). (4) 



Note that quantum discord is generally asymmetric un- 
der the change A-^r^ B, but in this paper we will consider 
only symmetric states for which the discord is naturally 
invariant under the choice of the measured subsystem. 

The hard step in the evaluation of the quantum dis- 
cord (4) is the optimization of the conditional entropy 
'H{ni}(^|-B) over all POVMs. Although this cannot be 
done in full generality analytically even for two qubits 
[26] the problem is tractable using analytical tools for 
certain subsets of states and POVMs, e.g., for Gaussian 
ones [27]. However, the question as to whether Gaussian 
measurements are globally optimal for the extraction of 
nonclassical correlations from general CV quantum states 
is essentially open and has not been settled even for Gaus- 
sian states themselves.^ 



B. CV Werner state 

We consider the CV Werner state [29] defined as 
p = p |V(A))(^(A)| + {l-p) p'}{ti) ® p'iifi), (5) 
where < p < 1, 



|V(A)) = Vl-A2^A"|n,n)^B. 



(6) 



n=0 



is the two-mode squeezed vacuum state with A = tanhr 
(r is the squeezing parameter) and 

oo 

is the thermal state with /x^ = {nj)/{l + {nj)), where {rij) 
is the mean number of thermal photons in mode j. We 
observe that the infinite-dimensional state (5) also pos- 
sesses the same structure as d-dimensional states invari- 
ant under the maximal commutative subgroup of U{d) 
introduced in [38]. 



III. SPECIAL CASE n^O 

First, let us consider the simplest special case of a 
Werner state with = which gives using Eq. (5) 

Po=p|V'(A))(^(A)| + (l-p)|00)(00|, (8) 



^ In this respect it is known that, for certain Gaussian states, 
non-Gaussian measurements are needed to minimize two-way 

quantifiers of nonclassical correlations such as the ameliorated 
measurement-induced disturbance [36]. For the quantum discord 
itself, Gaussian POVMS arc instead conjectured to be optimal 
among all possible CV measurements [27, 37]. 



representing just a mixture of a two-mode squeezed vac- 
uum state with the vacuum. For p > the partially 
transposed matrix p^^ (obtained by transposing po with 
respect to the degrees of freedom of subsystem A only) 
has negative eigenvalues [29] and therefore, according to 
the PPT criterion [34], the state (8) is entangled. Note 
that the state (8) has been further studied in [30, 39] from 
the point of view of its entanglement properties, as mea- 
sured by the negativity [40] , highlighting its applications 
for quantum key distribution. 



A. Exact calculation of quantum discord 

In order to calculate the entropies arising in the ex- 
pression of quantum discord (4) wc need to determine 
the eigenvalues of the reduced state po,B, the global state 
(8) and the conditional state pA\i = TTB[^B{i)Po]/PB{i)- 
The latter two states attain the form 



'^ = Cl|<^l)('/'l|+C2|</'2)(02|, 



(9) 



where Ci + Ca = 1 and |(/)i,2) are generally nonorthogo- 
nal normalized pure state vectors. The state (9) has at 

most two-dimensional support spanned by vectors \4'i^2) 
corresponding to eigenvalues 1^1,2 that read as 



l±v/l-4ClC2(l-|(</'l|^2)|^) 
^1,2 = ^ • (10) 

On inserting the eigenvalues (10) into the formula for the 
von Neumann entropy 



S{a) = - z/j In Vj 



(11) 



i=l 



we get analytically the entropy of the state (9). 

Returning back to the state (8) we get, in particular, 

|0i) = m\)),\4>2) = |00),Ci = p and C2 = 1 - p; the 
eigenvalues thus amount to 



z^i,2 = 



l±^/l-Ap{l-p)X^ 



(12) 



Hence, we can immediately calculate the entropy S{po) 
using formula (11), 



l-y/l-4p{l-p)\A / i-^i-4p(l-p)A^ 



(13) 

Tracing, further, the state (8) over mode A yields a di- 
agonal reduced state 



Po,s =PP^^(A) + (1-p)|0)b(0|, 



(14) 



with pg given in Eq. (7), possessing the eigenvalues 
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FIG. 1: Quantum discord V [Eq. (17)] versus the probabil- 
ity p and the squeezing factor A for the CV Werner state po 
[Eq. (8)]. AU the quantities plotted are dimensionless. 



2\ \2n 



i>o = l-p\ , i'n>o = P(l - A )A 



(15) 



for n e N. Making use of the definition (11) we get also 
the entropy of the reduced state in the form: 



S{po,b) = -|ln(l-pA2)+pA2ln 
2pA2 In A 



P(l-A^) 
1 - pA2 



+- 



1- A2 



(16) 



We are left with the minimization of the conditional 

entropy H{n,}{A\B) over all POVMs. The state (8) is 
constructed in such a way that we can "guess" the opti- 
mal measurement before attempting to solve this anyway 
hopeless task with brute-force. A close look at the state 
reveals that the optimal measurement on mode B is just 
photon counting (i.e., in this case, projection onto the 
eigenstates of po^s) characterized by the set of projec- 
tors {lis (to) = |to)b(to|}, where \m) is the m-th Fock 
state. Clearly, if we detect to photons in mode B in 
the state (8), the conditional state PA\m is simply a pure 
Fock state |to)^(to| which implies, immediately, that the 
conditional entropy attains its minimal possible value 
^{n(rTi)}(^l-S) = 0. This proves rigorously that photon 
counting is the globally optimal measurement strategy 
for the state po, thus giving the discord exactly equal to 



V{po)=S{po,b)-S{po), 



(17) 



where the involved entropies are defined in equations (13) 
and (16). The discord is an increasing function of both 
A and p and is plotted in Fig. 1. Note, that the consid- 
ered Werner state (8) belongs to the class of maximally 
correlated states for which Eq. (17) can be proved alter- 
natively [41] using the duality relation between classical 
correlations and entanglement of formation [42]. 



B. Discord versus Gaussian discord and the 
non-Gaussianity of the state 

A relevant question is whether photon counting is 
the only measurement realizing the global minimum in 
the evaluation of quantum discord, or, more specifically, 
whether there can exist Gaussian POVMs that attain in- 
stead an equally optimal performance on the state (8). 
The answer, as we argue in the following, is that nonclas- 
sical correlations in the state (8) captured by discord (17) 
cannot be extracted equally well by any Gaussian mea- 
surement. In other words, Gaussian discord (= 2?'^(/3o)), 
defined by Eq. (4) by imposing the minimization be re- 
stricted only to Gaussian measurements [27], is just a 
strictly larger upper bound on the true discord (17) for 
the mixed non-Gaussian CV Werner state (8). To show 
this, consider the following Gaussian POVM [43] consist- 
ing of elements 



where 



with ^ 



\a,£,) = D{a)S{m, 



(18) 



(19) 



tc^"^^ , is a pure normalized momentum-squeezed 
vacuum state with squeezing parameter t G [0,oo), 
that is rotated counterclockwise by a phase (f G [0, tt) 
and that is subsequently displaced by a G C. Here, 
D{a) = exp(aa^ — a* a) is the displacement operator and 
^{0 = { ~ /2} is the squeezing opera- 
tor. Note that t = Q corresponds to heterodyne detec- 
tion, whereas homodyne detection is obtained in the limit 
t oo. If the POVM clement 11(0:) is detected on mode 
B of the two- mode squeezed vacuum state (6), then mode 
A collapses into the state ]^,a; = se"'^^)^ which is a pure 
momentum-squeezed state with squeezing parameter 



1 



In 



1 



cosh(2r) 



cosh(2r) -H e^ 



(20) 



that is rotated clockwise by phase and that is displaced 

by 

P = [(-+ + z.)a* + (.+ - ._) e-'^^a] , (21) 

where z± = [cosh(2r) -I- exp (±2t)]~^. Thus, the ob- 
tained conditional state pA\a = Tr^ [lis {a) po] /q{a) is 
again a convex mixture of the form (9) where j^i) = 
1/3,0;), |</>2) = |0) and 



Ci 



pu{a) 
7rg(a) 



, C2 



{l-p)v{a) 
'Kq{a) 



(22) 



Here q{a) = \pu{a) + {1 — p)v{a)] /it is the probabil- 
ity density of obtaining the measurement outcome a, 
u{a) = (a, ^|/9^(A)|a, where Pb(A) is given in Eq. (7), 
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FIG. 2: Logarithmic plot of quantum discord (solid curve) and 
Gaussian quantum discord (dashed curve) versus the proba- 
bility p for the CV Werner state (8) with (from bottom to top) 
A = 0.1,0.5,0.9. All the quantities plotted are dimensionless. 



and v{a) = \{0\a,^)\'^. The latter function v{a) can be 
computed straightforwardly using the formula [44] 



■\/ cosh(t) 



(23) 



The overlap |(<;ii|(/'2)P = |(0|/3,w = se~'^'^)p appearing 
in eigenvalues (10) can be calculated exactly along the 
same lines. Expressing finally the thermal state Ps'(A) 
as 



1-A^ 
7rA2 



X'-' 



C 



10 (24) 



we get the function u{a) in the form 
1-A2 



u{a) = 



7rA2 



a-<;,0\'d'<;, (25) 



where we used the property of displacement operators 
D{-<i)D{a) = exp[(<^*a - ^a*)/2]D{a - c,). Using once 
again the formula (23) to express the overlap |(0|a— £) p 
and performing the integration over we arrive at the 
formula 



u{a) = 



l-A^ 



cosh(f)'^ 



X exp • 



1 - A4tanh^(t) 
(1-A2) [1 - Ananh2(t)] 



(26) 



A^ tanh^(t) 



(1 - A2)^tanh(t) 
2 [1 - A4tanh2(t)] 



(e-»2^a2+e^2^a*2^ 



Substituting the obtained explicit expressions for func- 
tions u{a)^v{a) and q{a) into Eqs. (22) and using the 
explicit expression for the overlap |((/)i|(/>2)P = |(0|/3,a; = 
gg-t2¥>^|2 from Eq. (10) the eigenvalues and hence 

the entropy S {^PA\a) of the conditional state pA\a- Sub- 
sequent averaging of the entropy over the density q{a) 



finally yields the Gaussian conditional entropy 

^fn(.)} (AB) = ^ q{a)S {pa\o) d'a (27) 

as a function of the squeezing parameter t and phase (f of 
the Gaussian measurement (18). Due to the complicated 
dependence of the conditional entropy S {pA\a) on a, the 
remaining integration over the complex plane C, where 
(Pa = d(Rea)d{lTna), has to be performed numerically. 
Likewise, minimization of the entropy (27) with respect 
to variables t and ip also requires numerics. This analysis 
reveals that, within the Gaussian POVM set, the entropy 
is minimized by homodyne detection on mode B. The re- 
sulting plots of Gaussian discord and the true quantum 
discord (17) (the latter obtained by photon counting on 
B) are shown in Fig. 2. The figure clearly shows that 
apart from trivial cases p = 0,1 the Gaussian discord is 
always strictly larger than the discord (17), meaning that 
general Gaussian measurements are strictly suboptimal 
(or, in other words, non-minimally disturbing) for the ex- 
traction of nonclassical correlations in the non-Gaussian 
state (8). 

To shine light on the reasons behind this, we compare 
the gap 



(28) 



quantifying how much Gaussian measurements overesti- 
mate the actual amount of nonclassical correlations, with 
a measure of the non-Gaussianity 6o of the state (8) [33]. 
The latter is defined as the quantum relative entropy be- 
tween po and a reference Gaussian state tq with the same 
first and second moments. 



So=S {to) - S (po) 



(29) 



where S (po) is given in (13). To find <S(to), we note 
that po has zero first moments but possesses a covariance 
matrix of the form 



To 



/C 5 \ 
C -5 
S C 

\0 -S c J 



(30) 



with C = pcosh(2r) + {1 — p), S = j)sinh(2r). The von 
Neumann entropy of the Gaussian state tq can be then 
written as [45] 



S{to) = + 



In 



(31) 



where v = v^[l - (1 - 2p)'^X^]/{l - X^) is the doubly- 
degenerate symplectic eigenvalue [46] of the covariance 
matrix Fg . The non-Gaussianity is a concave function 
of p and increases with A, diverging in the limit of infinite 
squeezing. In the regime of low squeezing, A ^ 1, a 
series expansion (up to the quadratic term in A) returns 
an approximate expression for the non-Gaussianity, 



c(A<l) 



(- 1 + p)pA2 [- 1 + ln(p( 1 - p) ) + 2 In A] . 
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(a) (b) 

FIG. 3: Normalized gap AOnorm [Eq. (33)] between Gaussian discord and optimal discord (dashed curve) and non-Gaussianity 
do (solid curve) of the CV Werner state (8), plotted as functions of (a) the probability p for different values of A (A = 0.2 and 
0.8 from bottom to top), and of (b) the squeezing factor A (at p = 0.5). The plots are in logarithmic scale. All the quantities 
plotted are dimensionless. 



On the other hand, the gap in discord just coincides 
with the conditional entropy (27) corresponding to ho- 
modync detection, AV = 'H'^Yi{ay}i^\-^)- low 
squeezing regime, we can also expand in series (up to 
the quadratic term in A) the integrand in (27), so as to 
obtain an approximate analytic expression for the gap, 

Ap(^«i) « 7r-^(-l+p)pA2[l-7-ln2+ln(p(l-p))+21nA], 

where 7 w 0.577 is Euler's constant. Defining the ratio 



ln(A)+7-l 



(32) 



between the approximate expressions for and AV (at 
p = 0.5), wc see that the linearly dependent relationship 
do « 4>aAP 

holds with good approximation for small A. In other 
words, for low squeezing (say A < 0.2), the (normalized) 

gap 



AV„ 



^xAV 



(33) 



between optimal Gaussian (homodyne) and globally op- 
timal non-Gaussian (photon counting) measurements 
for the extraction of nonclassical correlations, correctly 
characterizes and quantitatively reproduces the non- 
Gaussianity ^0 of the considered state (8). Interestingly, 

lim — — = <Po = IT . 
x^o AV 

This intriguing connection between the nonclassicality 
gap and non-Gaussianity fails to hold for larger values of 
A; the discrepancy between the two parameters becomes 
extreme in the limit A — > 1, when the non-Gaussianity 
diverges while the gap AV closes to zero. A comprehen- 
sive comparison between the normalized gap in discord 
(33) and the non-Gaussianity measure (29) is shown in 
Fig. 3. 



Finite versus infinite-dimensional POVMs 



One may argue, that even a simpler non-Gaussian state 
than that of given in Eq. (8) can be found possessing a 
strictly lower discord for a non-Gaussian measurement 
than for the best Gaussian measurement. For instance, 
the optimal measurement minimizing the discord in the 
qubit Werner state [28], studied in the seminal paper 
on quantum discord [16], is a simple non-Gaussian pro- 
jection onto the first two Fock states )0) and ]1). One 
can easily check that the optimization over all Gaus- 
sian measurements gives a strictly higher discord. Let 
us stress that the nonclassical correlations captured by 
discord are fundamentally different for the CV Werner 
state (8) considered here and for the qubit Werner state. 
Namely, although oiu- CV Werner [29] is formally a qubit- 
like state, the globally optimal POVM has an infinite 
number of elements given by projectors onto all Fock 
states. Moreover, it is not difficult to show that no 
POVM measurement on mode B possessing a finite num- 
ber N of elements fti, i = 1, . . . ,N can nullify the con- 
ditional entropy H{A\B) and hence also be globally op- 
timal. Namely, the conditional state corresponding to 
detection of the element fi^ has to be a pure state, that 
PA\i = Trs [po^i] = |x»)(Xi| in order for the entropy 
of the conditional state to vanish. Now, consider the 
element Qq = — X^ili ^i- The corresponding con- 
ditional state = TiB [po^q] = Pa.A - Ejli IXi)(Xi|, 
where po.A is obtained from Eq. (14) by replacing B with 
A, cannot have neither zero, nor one, nor even any fi- 
nite number of strictly positive eigenvalues, as this would 
imply that the state po.A also has a finite number of 
strictly positive eigenvalues which is not the case [see 
equation (15)]. Therefore, for any finite N the condi- 
tional state p^iQ is definitely a mixed state possessing 
strictly positive entropy and hence resulting in a strictly 
positive and therefore suboptimal conditional entropy 
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(i.e., non-optimal discord). Consequently, the globally 
optimal POVM of the qubit Werner state is only two- 
component and thus the qubit Werner state is only a 
trivial embedding of a two-qubit state into an infinitely- 
dimensional two-mode state space carrying only qubit- 
type non-Gaussian nonclassical correlations. In con- 
trast, the CV Werner state (8) carries genuinely CV non- 
Gaussian nonclassical correlations that can be optimally 
extracted only by a non-Gaussian POVM measurement 
with an infinite number of elements: in this particular 
case, photon counting. 



D. Extension to mixtures of n Gaussian states 

Before going further let us note that, for a state of the 
form (9), the global optimality of photon counting (for 
the calculation of quantum discord) follows immediately 
from the fact that the projection of one of its modes onto 
a Fock state projects the other mode onto the same Fock 
state. As a consequence the conditional entropy achieves 
the minimum possible value 7i{u{m)}{A\B) = and the 
discord is of the form (17). It is not difficult to find 
more general mixed states with the same property, for 
instance, states with the structure 

oo 

Pg = X] 1mn\mm){nn\, (34) 

m,n=0 

with g*„ = q-mn, J2m=Q 1-mrn = 1 and the matrix Q with 
elements Qij^ki = Qik^ij^ki being positive-semidefinite. 

A particular example of such states is a convex mix- 
ture of an arbitrary number of two-mode squeezed vacua 
(6) with different squeezing parameters obtained for 
Qmn = J2iP'i (l ^ ^i) where pi are probabilities 

and Aj Xj for i ^ j. We remark that for all non- 
Gaussian states of the form (34) the quantum discord 
can be computed exactly. 



E. Other noncleissicality indicators 

The quantum state (8) admits furthermore an analyt- 
ical calculation of other optimized entropic quantifiers 
of nonclassical correlations [19], encompassing the 'ame- 
liorated' measurement- induced disturbance [22] and the 
relative entropy of quantumness [24, 25]. Interestingly, as 
we will now show, both of them coincide with the discord 
(17). This is in complete analogy with the correspond- 
ing Werner state for two qubits [20, 22, 26] and for two 
qudits with arbitrary d [47]. 



1. Measurement-induced disturbance 

The optimized nonclassicality indicator called 'ame- 
liorated' measurement-induced disturbance (AMID) is 



a two-way measure of quantum correlations defined as 

[21, 22, 36] 

■^{Po) ^Iq{po) -^c{pa), (35) 

where I, (po) = >5(po,^)+>5(po,b)->5(po) is the quantum 
mutual information and 

Ic{po) = sup I{pab) (36) 

is the classical mutual information [48] of the quantum 
state po- Here I{pab) = 'H{pa)+'H{j?b)-'H{pab) is the 
classical mutual information of the joint probability 6Ss- 
tribution pAB(fc, = Tr[poIlA(fc)®nB(Z)] of outcomes of 
local measurements 11^ and lis on where ^{pab) and 
'H{Pa) {H{pb)) are the Shannon entropies of the joint dis- 
tribution and marginal distribution p^i (fc) = '^iPAB{k,l) 
{pb{1) = J2kPAB{k,l)), respectively. The classical mu- 
tual information is upper bounded as [21] 

I{pab) < mm{S{po,A),S{po,B),Iq (po)} ■ (37) 

Since S{po^a) = S{po,b) for the state (8), we need to 
compare 5(po,s) with Iq (po) which can be done using 
majorization theory for infinite-dimensional density ma- 
trices [49] . Consider two such density matrices A and B 
with ai, 02, . . . and 6i, 62, . • . being their non-zero eigen- 
values arranged in a decreasing order and repeated ac- 
cording to their multiplicity. If some of the matrices, for 
example, A has only a finite number k of non-zero eigen- 
values, we set ttfe+i = afe-i-2 = . . . = 0. We then say that 
A is more mixed than B, and write Ay B, ii 

k k 

Y,ai<Y,b^, fc = l,2,.... (38) 

i=l i=l 

It holds further [49], that ii A y B then their von Neu- 
mann entropies satisfy S{A) > S{B). Taking now the 
eigenvalues Dj, j = 0, 1, . . . given in Eq. (15) instead of 
eigenvalues 04, and 1/1,2 given in Eq. (12) instead of eigen- 
values bi, one easily finds that they satisfy Eq. (38). This 
implies that the density matrices po and po,B of Eqs. (8) 
and (14) satisfy po,s >- Po and therefore »S(po,b) > ^{Po)- 
Hence, one gets Xg (po) > 5(po,b) which leads, using the 
inequality (37), to the upper bound on the classical mu- 
tual information in the form X{pab) < S{Po,b)- Consid- 
ering now photon counting on both modes in the state 
(8), one gets 'H(pa) = 'H(pb) = H(pab) = 5(po,a) = 
'5(/Oo,-b) which gives the classical mutual information in 
the form I{pab) = ^(po.b)- Hence, the latter inequal- 
ity is saturated by photon counting which finally yields 
A{po) = 5(po,b) - S{po) = V{po), that is, AMID coin- 
cides with the quantum discord, Eq. (17). 

2. Relative entropy of quantumness 

The measure quantifies the minimum distance, in 
terms of relative entropy, between a quantum state po 
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and the set of eompletely classically correlated states [24] . 
It is operationally associated to the amount of distillable 
entanglement that can be generated between a quantum 
state po (carrying nonclassical correlations) and a set of 
ancillary systems in the worst case scenario of an activa- 
tion protocol [25]. The relative entropy of quantumness 
is defined as [24, 25] 

Q (po) = min ^ [H {pab) - S (po)] , (39) 

where the minimization is performed over local mea- 
surements Ha and Hb consisting of collections of one- 
dimensional projectors {n^(fc)} and {TIb{1)}- 

The minimum in Eq. (39) can be found using the same 
trick as in the previous cases. First, we take a suit- 
able tight bound on the quantity that is to be optimized 
and then we guess a measurement saturating the bound. 
In the present case we look for the minimum of the 
joint Shannon entropy V. {pab) satisfying the inequalities 
H {pab) > "H {pb) > S {pq,b)- Using the previous result 
that for photon counting we get T-L{pab) = S{po^b), we 
see that the measurement saturates the lower bound and 
therefore Q (po) = S{po.b) — S (po). Thus all three con- 
sidered nonclassicality indicators coincide for the state 
(8), i.e., 

V{po) = A (po) = Q (po) = S{po,b) - S{po). (40) 



IV. GENERAL CASE 

Let us now move to the analysis of nonclassical cor- 
relations in a generic CV Werner state (5) with p 7^ 
0, complementing the seminal analysis of nonclassical 
correlations in a two-qubit Werner state performed in 
[16, 20, 22, 26]. This can be interesting in particular be- 
cause, in contrast to the qubit case, there can potentially 
exist PPT entangled CV Werner states [29]. 

In the present general case we do not have any tight 
bounds, similar to those of the previous Section, allow- 
ing us to perform exact optimizations in Eqs. (4), (36) 
and (39). For this reason, we cannot prove the global 
optimality of photon counting or any other measurement 
strategy analytically. We then resort to computing up- 
per bounds on discord, AMID'^ and relative entropy of 
quantumness, obtained for (possibly nonoptimized) mea- 
surements in the local eigcnbasis of the reduced state (s) 
of the two-mode CV Werner state. Interestingly, all the 
upper bounds on the different quantities again coincide 
as we show later in this Section: this hints at the conjec- 
ture that they might be indeed tight for the considered 
states, although we cannot provide conclusive evidence 
of this claim. We also derive nontrivial lower bounds for 
the nonclassical correlations. 



A. Upper and lower bounds on discord 

We consider a nonoptimized upper bound on quantum 
discord defined for a density matrix p as 



U{p)=S{pB)-S{p)+n,i^{A\B), 



(41) 



where 'Hcig{A\B) is the conditional entropy for the mea- 
surement of mode B in the local eigenbasis of the reduced 
state pB (see also [20]). For the general CV Werner state 
of (5), tracing p over mode A gives the reduced state 



PB=pp'^{X) + {l-p)p'B^{^i), 



th 



(42) 



where p"^ is defined in Eq. (7). This state is diagonal in 
the Fock basis with eigenvalues 



gm=p{l- X') A^'" + (1 - p) (1 - p2 



(43) 



that give, after substitution into Eq. (11), the marginal 
entropy S{pb) appearing in Eq. (41). 

The local eigenbasis is a Fock basis and so the projec- 
tion on it corresponds again, even in the present general 
case, to photon counting. The conditional state PA\m — 
Trs [lTO)B(mlp] /pB(m), where pB{m) = B{m\pB\m)B, 
obtained by projecting mode B onto Fock state ]m) reads 
explicitly 



PA\t 



p (1 - A2) \^'^\m)A{m\ + (1 -p)(l - p')p2'"p^^ (p) 



PB{m) 



(44) 



with pB{m) = p (1 - A2) A^"* -t- (1 - p)(l - p2)^2m_ 
has the eigenvalues 



'In 



2m 



5mn + (l-p)(l-p^)>' 



2\2 2{m+n) 



PB{m) 



(45) 



where Smn is the Kronecker symbol, that give the follow- 
ing entropy of the conditional state 



>5 {pA\m) = - 



(l-p)(l-p2)V2 



X <^ In 



Pb (m) 

(1-p) (l-p2)^ 



PB{m) 

(^-p-)+ln(p^) 

_l-p2 (1-^2)2 
7m im 



(46) 



^ In this case the upper bound on AMID is simply the nonopti- 
mized measurement-induced disturbance (MID) [20]. 



Hence, one gets the conditional entropy 'Heig(^l-B) = 

{pA\m)- 

It remains to calculate the global entropy of the state 
(5). For this purpose it is convenient to express the state 
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as 



/9= ^ Mmn\ni,m){n,n\+ ^ e„„|m,n)(m,n| 

m^n=0 



m,n=0 



where 



(47) 



+ (48) 

emu = (1-p)(1-m')%'^'"+"^- (49) 

The state (5) thus possesses the eigenvalues Cmn corre- 
sponding to the eigenvectors |m,n), m ^ n = 0,1,... 
and the remaining eigenvalues (= /;) are the eigenvalues 
of the infinite-dimensional matrix M with elements (48). 
This gives the global entropy 

^(p) = -^4fe^{in[(i-P)(i-.^r] 



+ 



21n(M)(l + M^ + 2/i^) 



(50) 

The eigenvalues /; of matrix M appearing in the last ex- 
pression of the previous equation cannot be calculated 
analytically and one has to resort to numerical diagonal- 
ization of a sufficiently large truncated matrix. Hence, 
one gets using Eq. (50), and expressions for T-Leig{A\B) 
and S{pb), the sought upper bound (41) on the true 
quantum discord. 

The true discord can be also bounded from below 
in the following way. Let us observe first, that apart 
from the trivial case p = (corresponding to a prod- 
uct state) all other CV Werner states have nonclassical 
correlations as they possess a strictly positive quantum 
discord I'(p) > 0. This can be proven using the sufR- 
cient condition on strict positivity of quantum discord 
[50] according to which Dip) > for a state p if at least 

one off-diagonal block pl- = B{i\p\j)B, * 7^ J is not 
normal, i.e., it does not commute with its adjoint. In 
the present case of the Werner state (5) we have explic- 
itly p^f^ = p (1 — A^) A*+-' |z)B(j|. Assuming p > and 
< A < 1 this gives immediately a nonzero commutator 



^0 for i ^ j, p ^ 0, 



\j)b{j\] 

(51) 



as required. 

An explicit, non-tight lower bound can be derived that 

is nonnegative (and thus nontrivial) at least on some 
subinterval of probabilities p. Namely, assume a POVM 
on mode B given by a collection of rank-1 operators 
{\'4>j){'4'j\}- If the component \'4>i){'4'i\ is detected on 
mode B in the state (5), then mode A collapses into the 
normalized conditional state 

PA\i = , (52) 

Pi 



where \4>i) a is a pure unnormalized state that is not spec- 
ified here and pi = {ipi\pB\'4'i) with pB given in Eq. (42) 
is the probability of measuring the outcome i. The state 
is a convex mixture of a pure state and a thermal state. 
Making use of the concavity of the von Neumann en- 
tropy S (^jPjPj^ > J2jPj^ (Pj) ^^'^ the fact that the 
entropy vanishes on pure states we arrive at the following 
inequality: 



•5 {pA\i) > 



Pi 



(V;iK(M)IV'i)5[p^^M)]- (53) 



By multiplying both sides of the inequality by pi and 
summing over i one finds the classical conditional en- 
tropy to be lower bounded as H{\^.^^.\-j{A\B) > (1 — 
p)S [p*^ (p)] which yields finally using Eq. (4) the lower 
bound 

£(p) = S (pb) -S{p) + {l-p)S [p'^ ip)] . (54) 

In what follows we evaluate the upper and lower 

bounds in Eqs. (41) and (54), respectively, for two par- 
ticularly important two-parametric subfamilies of the set 
of CV Werner states. 



B. Examples 

1. The case X = p 

First, we consider the case X = p. Then the reduced 
state (42) is just a thermal state pb = Pb (A) with a 
well-known entropy 



(1-A2) 1-A2 \1-AV 



(55) 



In this case the CV Werner state is a mixture of a pure 
two- mode squeezed vacuum state (6) with a product of its 
marginals; in the strong squeezing limit A — > 1, this state 
approaches a mixture of a maximally entangled EPR 
state and a maximally mixed (infinitely thermal) state, 
which is a direct CV counterpart to the qubit Werner 
state [28]. Further, it can be shown that the CV Werner 
state (5) with A = /i > is entangled for any p > [29]. 
The upper bound (41) and lower bound (54) on quan- 
tum discord are depicted in Fig. 4. Note that in this and 
in the following plots, only nonzero values of the lower 
bound (54) will be shown. 



2. The case X = p* 

The case for which A = /it^ is interesting because the 

CV Werner state (5) runs through three different sepa- 
rability regions as the parameter p increases [29] : 

1. If p < Pscp = 2(i-^24_^4) , then the state p is sepa- 



rable (dark gray strip in Fig. 5). 
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2.0 




FIG. 4: Upper bound U, Eq. (41) (solid curve), and lower 
bound C, Eq. (54) (dashed curve), on quantum discord (4) 
versus probability p for the CV Werner state (5) with A = 
/i = 0.8. All the quantities plotted are dimensionless. 

2. If pscp < p < pppT = (i_^2)2^(ii^8)^2 . the state 
p is PPT (i.e., it has positive partial transposition 
[34]) and it is unknown whether it is undistillablc- 
entangled or separable (light gray strip in Fig. 5). 

3. lip> pppt, then the state p is non-PPT and there- 
fore entangled (white region in Fig. 5). 

C. Other nonoptimized nonclassicality indicators 

Next we focus on the determination of the nonop- 
timized (upper bound) version of AMID (35) called 
measurement-induced disturbance (MID) defined as [20] 

M{p)=Ig{p)-X{pAB), (56) 

where I{pab) is the classical Shannon mutual informa- 
tion of a probability distribution of results of measure- 
ments in eigenbases of the reduced states Pa,B: which in 
the present case coincides with the joint photon-number 
distribution pyiB(m, n) ~ AB{mn\p\mn) ab- Its marginal 
distributions pa = Pb coincide with the eigenvalues of 
the reduced states (43), i.e., PA{m) = PB{m) = Qm 
whence we get the equality between local Shannon and 
von Neumann entropies 

H{pa)=H{pb)=S{pa)=S{pb). (57) 
Hence MID simplifies to 

M{p)=n{pAB)-S{p). (58) 

The global Shannon entropy can be derived easily by 
noting that the eigenvalues (45) of the conditional state 
(44) satisfy r^l™'' = pAB{m,n)/pB{m) thus representing 
a conditional probability PAB{n\m) of detecting n pho- 
tons in mode A given m photons have been detected in 
mode B. This implies immediately that 'Hcig(^l-B) = 
{Pab) — S{pb), where we have used Eq. (57) and the 




FIG. 5: Upper bound U, Eq. (41) (solid curve), and lower 
bound C, Eq. (54) (dashed curve), on quantum discord (4) 
versus probability p for the CV Werner state (5) with \ = p,* 
and p = 0.8. The dark gray shaded region corresponds to 
separable states, the light gray shaded region corresponds to 
PPT states with unknown separability properties, and the 
white region corresponds to entangled non-PPT states. The 
boundary probabilities psep and pppT are depicted by vertical 
solid line and dashed line, respectively. All the quantities 
plotted are dimensionless. 



equality pA{m) = pB{m). Substituting from here for 
{Pab) into Eq. (58) finally leads to the equality of the 
upper bound on discord (41) and MID (56) 

M{p)=U{p). (59) 

Note that the two coincident quantities also provide 

an upper bound for the relative entropy of quantum- 
ness [19, 24] of the states (5). Note further that the 
lower bound (54) on discord is also a lower bound for 
the other measures of nonclassical correlations such as 
AMID, MID, and relative entropy of quantumness, since 
quantum discord is in general smaller than those men- 
tioned quantities for arbitrary bipartite quantum states 
[19, 21]. 

V. PARTIALLY TRANSPOSED CV WERNER 
STATE 

One of the main technical disadvantages of the CV 

Werner state with /i 7^ is that its eigenvalues, and con- 
sequently its global von Neumann entropy, cannot be cal- 
culated analytically. Interestingly, this ceases to be the 
case if the state is partially transposed [29]. Then, one 
can find regions of parameters p, X and p for which the 
partial transposes are positive-semidefinite, and so repre- 
sent a legitimate quantum state in their own right. Thus, 
one obtains another, sometimes simpler to treat family of 
non-Gaussian quantum states for which we can get fur- 
ther with analytical tools than in the case of the original 
state. 

Let us illustrate this on a simple example of the CV 
Werner state (5) with X = p^. As was shown in Ref. [29] 
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for 



P 



1 A 



(60) 



the partial transposition p^^ = p oi the Werner state p 
with respect to mode A, 

p = M A'"+" (|n, m)(m,n| + |m, n)(m, n|) , 

m,n— 

(61) 

possesses the following nonnegative nonzero eigenvalues 



2NX 



m-\-n 



m = 0,l,..., (62) 
m > n = 0, 1, . . . , (63) 



where A/" = (l — A^) (1 — A) /2, and is therefore a valid 
two-mode density matrix corresponding to a different 
non-Gaussian state. Direct substitution of the eigenval- 
ues into Eq. (11) gives the analytical expression for the 
global entropy of p of the form 



S{p) = - 



1 I Q \ 



(64) 



Tracing the state (61) over mode A one gets the reduced 
state of mode B 



m=0 



A" 



PB=MY. (a2- + — )|m)B(m| 
with entropy 



(65) 



S{~Pb) 



AA^ A 



TO=0 



1 - A 



In A 



1 



, A(1 + 3A), / 

+ ^^(^)+2(l^^^^. 



1- A, 
(66) 



Similarly one can find a reduced state pA of mode A 
which coincides with the reduced state (65) and yields 
the local entropy S {pa) = S {pb)- 

Upon detecting m photons in mode B in the state (61), 
mode A collapses into the normalized conditional state 



PA\m 

where 



Psim) 



A2'"|m)A(m|+A'"^A'=|A:)A(fc|j , 

(67) 



fe=0 



ps(m) =A/'A'" A'" + 



1-A 



(68) 



is the probability of detecting m photons on mode 

B in the state (61). After some algebra, the corre- 
sponding nonoptimizcd conditional entropy 'Heig(^l-B) = 
^^^qPb(to)'5 {pA\rn) attains then the form 



Heis{A\B) = S (p) - S (pb) + Xln2, 



(69) 



0.7 r 

0.6^ 

0.5; 

S 0.4 ; 

0.3: 

0.2 
0.1 
0.0^ 



0.0 0.2 0.4 0.6 0.8 1.0 

A 

FIG. 6: Upper bound U, Eq. (70) (solid curve), and lower 
bound £, Eq. (71) (dashed curve), on quantum discord (4) 
versus squeezing parameter A = tanh r for the partially trans- 
posed CV Werner state (61). All the quantities plotted are 
dimensionless. 



where we have used Eqs. (64) and (66). Substituting 
finally the latter formula into Eq. (41) we arrive at a 

very simple analytical expression for the upper bound on 
the quantum discord of the state (61), 



U{p) = Aln2. 



(70) 



Even for the partially transposed CV Werner state p, it 
is possible to derive a nontrivial lower bound on quantum 
discord (4). Repeating the algorithm leading to Eq. (54) 
for the state (61), one gets the lower bound in the form: 



C{p)=S{pB)-Sip) 



1 + A 
2 



S 



Pa 



(71) 



The upper bound (70) and lower bound (71) on discord 
are depicted in Fig. 6 as a function of the parameter A. 
Note that in this case they are quite close to each other, 
with the lower bound being always faithful in the whole 
considered parameter space. Also, nonclassical correla- 
tions in this PPT state (which may be separable or at 
most contain undistillable entanglement) are quite weak 
(yet always nonzero), increasing slowly with the squeez- 
ing r (recall that A = tanh r) and converging to the small, 
finite value In 2 in the limit r — )■ oo. 

As for other nonclassicality indicators, we get equiva- 
lent results. Moving for instance to the evaluation of the 
MID (56), one gets the joint photon-number distribution 
for the state (61) to be 

PAB (m, n) = ArA™+" (1 + S^n) ■ (72) 

The Shannon entropy of the distribution reads 

n (pab) ^ S ip) + Xln2, (73) 

where the global entropy <S {p) is given in Eq. (64) . The 
marginal distributions on each mode coincide and they 

are given by Eq. (68). Hence for the state p also, the local 
Shannon and von Neumann entropies satisfy Eq. (57). 
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Making use of the latter equality in the definition of MID 

(56) we can express it as 

M{~p)=H{vAB)-S{p)=\\n2, (74) 

where in the derivation of the second equality we used 
Eq. (73). Thus, as for the generic CV Werner state of 
the previous Section, for the considered partially trans- 
posed CV Werner state the MID coincides with the upper 
bound on discord associated with local photon counting, 
i.e.,M{~p)=U{~p). 

VI. CONCLUSIONS 

In this paper we studied the nonclassicality of correla- 
tions in a class of non-Gaussian states of a two-mode CV 
system. We adopted primarily quantum discord [16, 17], 
as well as other quantifiers such as measurement-induced 
disturbance and relative entropy of quantumness, and 
provided exact results (when possible) and in general up- 
per and lower bounds on their quantification in the con- 
sidered states. Our analysis enabled us to venture beyond 
Gaussian states and operations to access general features 
of quantum correlations in infinite-dimensional states. 
We focused in particular on two-mode CV Werner states 
[29] , constructed as mixtures of Gaussian states and thus 
possessing a positive, yet non-Gaussian-shaped Wigner 
function on the quantum phase space. Such states pro- 
vide important testbeds for the understanding of non- 
classicality beyond entanglement and its interplay with 
entanglement itself. For the special case of a two-mode 
squeezed state mixed with a vacuum [30], we were able 
to calculate quantum discord exactly, proving that local 
photon counting constitutes the globally optimal strat- 
egy to extract nonclassical correlations, and any Gaus- 
sian measurement strategy turns out to be suboptimal 
for the task. The considered states constitute probably 
the simplest example of bipartite states possessing gen- 
uinely non-Gaussian CV nonclassical correlations beside 
entanglement. The relationship between nonclassicality 
and non-Gaussianity was further highlighted by observ- 
ing that the gap between the quantum discord and its 
(non-optimal) counterpart, restricted to Gaussian mea- 
surements only, scales linearly with an entropic measure 
of non-Gaussianity [33] for the considered class of states 
in the regime of low squeezing. For general CV Werner 
states, photon counting measurements provide in general 
upper bounds on quantum discord and related nonclassi- 
cality measures, but we could not prove the tightness of 
such bounds analytically. We are tempted to conjecture 
that our upper bound does yield the true quantum dis- 
cord for all CV Werner states, and future progress on vali- 
dating or disproving this claim would be valuable. On the 



other hand, we derived nontrivial lower bounds on quan- 
tum discord as well for the considered states, which we 
do not expect to be tight. We finally constructed a par- 
ticular instance of a non-Gaussian state which is positive 
under partial transposition, and whose upper and lower 
bounds on discord are analytically computable and allow 
us to pin down its nonclassical correlations as being quite 
weak, staying finite even in the limit of infinite squeezing. 
Our study evidences a trend for quantum correlations to 
be generally limited in the absence of (distillable) entan- 
glement, as originally noted for Gaussian states [27]. It 
would be interesting to provide tight upper bounds on 
the attainable amount of quantum discord for all sepa- 
rable CV states. Such a bound is known for two-mode 
Gaussian states and corresponds to one unit of discord, 
on a scale ranging to infinity [27]. On finite-dimensional 
two-qudit systems, on the other hand, it is known that 
separable mixed states can have as much nonclassical cor- 
relations as maximally entangled pure states, approach- 
ing the bound In d when measured by the relative entropy 
of quantumness [24]. 

Nonclassicality and non-Gaussianity are two of the 
most important resources for the optimization and real- 
istic implementation of present-day and next-generation 
quantum technology [10, 11]. This paper realizes a first 
step to explore the interplay between the two in physi- 
cally relevant CV states. An important next target for 
future work would be to study the structure of non- 
classical correlations in other practically useful states 
deviating from Gaussianity, such as photon-subtracted 
states [13, 39, 51]. Comparing the performance of non- 
Gaussian measurements, such as photon counting, with 
that of Gaussian strategies, such as homodyne and het- 
erodyne detection, for accessing nonclassical correlations, 
and studying how the gap between the two scales with the 
non-Gaussianity of the states [33] , and with other non- 
classical parameters widely used in quantum optics [52], 
could provide novel insight into the nature of quantum- 
ness (in its broadest sense) and its potential exploitation 
for CV quantum information. 
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